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Abstract 

Any positive-energy state of a free Dirac particle that is initially 
highly-localized, evolves in time by spreading at speeds close to the 
speed of light. This general phenomenon is explained by the fact 
that the Dirac evolution can be approximated arbitrarily closely by 
a quantum random walk, where the roles of coin and walker systems 
are naturally attributed to the spin and position degrees of freedom of 
the particle. Initially entangled and spatially localized spin-position 
states evolve with asymptotic two-horned distributions of the position 
probability, familiar from earlier studies of quantum walks. For the 
Dirac particle, the two horns travel apart at close to the speed of light. 
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1 Introduction 



The concept of a quantum random walk (QRW) has been widely discussed 
and extended in various directions [I] since its introduction El Hj and 
developement 00IZ1IHII3IIDIIII1II21- Much of the interest has derived 
from an expectation that such a mathematically attractive idea should have 
important applications in quantum information theory, analogous to known 
applications of classical random walks (CRWs) in classical information sci- 
ence. 

On the other hand, CRWs also have many applications outside classical 
information theory, in a wide variety of areas of science where mathematical 
modelling is involved, so it should not be surprising if QRWs find applications 
outside quantum information theory. Here we describe such an application 
to relativistic quantum mechanics; for an earlier application in this field, see 

The evolution in time of the state of a free Dirac particle, starting from 
a highly localized, positive-energy state, is a quantum process that has only 
recently been described fully ^3] • There had been a widespread misapprehen- 
sion that no relativistic particle with nonzero rest-mass m could be localized 
much within its Compton wavelength Ac = h/mc, where c is the speed of 
light. However it has now been shown ^2] that there is no such difficulty 
for the Dirac particle if localization is characterized in terms of the Dirac 
position operator x, by making A x = a/ (x 2 ) — (x) 2 small while keeping the 
energy positive, and not by unrealistic attempts to restrict the domain of the 
wavefunction to a bounded region in configuration space. Arbitrarily precise 
localization, with A x <C Xc, is possible in the case of the free Dirac particle 
with positive energy. When the particle is localized in such an initial state, it 
has an associated uncertainty in energy A E ^> mc 2 , and the subsequent evo- 
lution produces a probability density that spreads outwards in all directions 
at close to speed c. The graph of the evolving density along any axis through 
the centre of initial localization (see Fig. 1 in |14j ) shows a striking resem- 
blance to the two-horned density found for a typical 1-dimensional QRW [I] . 
For the Dirac particle, the horns are close to distance ct from the starting 
point. We shall show that this is not a coincidence, and that the evolution 
of any positive-energy state of a free Dirac particle moving in 1-dimension 
can be modelled arbitrarily closely as a QRW of the type described in detail 
by Ambainis et al. [\ , Konno ^01 CZ] and others. 

In addition to providing a somewhat surprising application of a QRW to 
a real process, this result provides some new insights as to the nature of the 
quantum walk itself. Until now the various proposed realization schemes for 
QRWs were typically based on the idea that the coin and walker degrees of 
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freedom of the walk should be associated with two distinct quantum sys- 
tems. These two systems were to be combined by means of some form of 
dynamical coupling-decoupling scheme. The present application shows that 
alternatively, a single quantum mechanical object such as the free Dirac par- 
ticle — by its very nature as a relativistic system with translational and spin 
degrees of freedom — can be identified in the course of its time evolution 
with a quantum random walk. This natural occurence of a QRW, instead 
of some engineered realization, suggests that the question of its ontological 
status is still an interesting and open one. 

The present application draws attention to two other important features 
of QRWs that have been emphasized by others [HI CHI- The first is that 
a QRW is a unitary evolution; the associated randomness is of the kind 
associated with every unitary quantum evolution. In particular, a QRW is 
typically reversible in time, unlike a CRW. The second feature is that a QRW 
is typically a ballistic process, associated with spreading at a constant speed, 
unlike diffusive CRWs, where spreading is proportional to the square root 
of the time. These two features are essential for the application that we 
describe below to the free Dirac evolution, which is a time-reversible process 
characterized by spreading near the speed of light. 

In what follows, we relate the evolution of Dirac's equation to that of 
a QRW based on the canonical Heisenberg algebra extended by the Dirac 
matrices. Then we construct and discuss the limiting probability distribution 
describing the translational spreading of an initial state. The paper concludes 
with some speculations about the physical reality of the quantum walk of the 
Dirac particle, and the possibility of detecting it experimentally. 

2 QRW and free Dirac evolution 

The free Dirac Hamiltonian operator for a particle with zero momentum 
along the y and z directions is 



acting on 4-component spinor wavefunctions \fr(x). Here we adopt a repre- 
sentation of the Dirac matrices with 



where crj, i — 1, 2, 3 are the usual Pauli matrices, and I2 is the 2x2 unit 
matrix. In this representation the helicity (spin) operator associated with 
rotations about the rc-axis is £ = 1 2 <8> 03. From this point onwards we adopt 



H(p) = cap + mc 2 (3 , p = —ihd/dx , 




a = cr 3 ® cr 3 , (3 = a 2 <g> I2 



(2) 
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the natural units K — c — m — 1. Recalling that only those solutions of 
Dirac's equation with positive energy describe physical electron states, we 
introduce the orthonormal positive-energy spinors in momentum space 

1 /1 + E(p)±p 

U±(P) ~ 2^E{p){E{p) + l] [ - (1 + E{p) Tp) 

where e + = (1,0) T , e_ = (0, 1) T and E(p) = ■sjp 1 + 1. These spinors satisfy 
the relations 




u±(p) ] u±{p) = 1 , u±{p) ] u T (p) = 



H{p)u±{p) = E(p)u±(p) , Eu±(p) = ±-u±{p) . (4) 

Now we can write an arbitrary positive-energy wavefunction (with zero y and 
z components of momentum) ^q{x) in terms of two arbitrary functions f±{p) 
as 

¥ o (a0 = "4= r e ipX {/+(P)«+(P) + /-(P)«-(P)} • (5) 



/2tt y_ 

Then 



*o(a:) + *o(ar)^ = l<— > / {\U(p)\ 2 + \f-(p)\ 2 } dp = 1 . (6) 

Suppose now that we choose a normalized positive-energy state with a 
definite helicity +1/2 and finite mean energy. Then 



/+(p) = /(p), /-(p) = 0, / |/(p)| 2 *=l, (T) 

and 

/oo 
E(p)|/(p)| 2 dp = E <oo. (8) 
-oo 

With f-(p) = 0, the action of H(p) in the second factor of the tensor 
product space in © and Q becomes trivial, as the second spinor remains 
constant at the value e + . Thus the second factor space can be ignored, and 
we can consider H(p) to have the effective form 

H(p) = a 3 p + a 2 (9) 
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acting in the first factor space. In this first space, we write 

(1 0f=|+) (0 1) T =R, (10) 
so that the positive-energy spinor u + (p) in (JHJ) takes the (effective) form 

= o / T-i / W L , =ff {(! + #(p) + P)l+) + *(! + - P)|-» • 

(11) 

Next we consider a fixed, small time interval At -C 1/Eq. The (effective) 
unitary evolution operator for the Dirac particle can then be approximated 
over the time interval At using the relations 



e -iH(p)At = VU + ([At] 2 ), 

V = e ~iAta 3 p^ U = e -iAta 2 _ ^ 

Here we see the appearance of the evolution operator V U for a 1-dimensional 
QRW [3], with V enacting a step of length At to the left or right along 
the x-axis (the "walker space"), depending on the sign of <r 3 , and with the 
reshuffling matrix U representing the "quantum coin toss" after each time 
interval of duration At. For a longer time t = nAt, we have from (fT2|) 

e- iH{p)t = {V U) n + 0{At) , (13) 

and we see that the evolution of the state of the Dirac particle over any finite 
time t can be obtained arbitrarily accurately by replacing the exact evolution 
operator by (V U) n and letting n — > oo and At — * with nAt = t. In other 
words, 



lim ( e -^P e -^A n = e -iH(p)t (M) 

At^O 
nAt = t 

and we can emulate the Dirac evolution by the evolution of a QRW. If we 
are interested in particular in the evolution of the probability density along 
the x-axis for the electron, we need only go at each time t to the asymptotic 
form of the QRW probability distribution for the "walker" EE EI • 

It is important to note at this point that whereas the exact Dirac evo- 
lution operator e~ tH ^ t obviously preserves the positive energy condition for 
physically meaningful states, the same is not true of the approximate, QRW 



5 



evolution (VU) n . However, (jl3j) and (jl4j) show that in the asymptotic limit 
described, the positive energy condition is respected. 

We close this section with the following remark. Rewriting the evolu- 
tion operators as V = |+)(+| e~ iAtp + |-)(-| e iAtp , U = e - iAta2 , we identify 
the type of quantum walk involved here as a Canonical Algebra QRW (CA- 
QRW) in the classification of [HI] • In contrast to the Euclidean QRW, which 
takes place on the integers and whose evolution operator is constructed from 
the generators of the Euclidean algebra jTHl EHl EU E2j , in the present case 
the generators of the canonical Heisenberg algebra — position and momen- 
tum operators - - are used in the construction of a discrete walk on the 
x-coordinate axis. The close algebraic relationship between these two walks 
facilitates the solution of the time evolution in the present case, provided (as 
is done in next section) that we carefully discretize the coordinate-space (gen- 
eralized) eigenf unctions, which unlike their Euclidean QRW counterparts, are 
not orthogonal. 

3 Asymptotic solutions and localization 

Let Tt denote the Hilbert space spanned by all vectors |$) <g) |±), correspond- 
ing in the coordinate representation to normalizable 2-component wavefunc- 
tions $(x)|±). Introduce a dense subspace S < Tt consisting of all finite 
linear combinations of suitably regular vectors |$) ® |±) G Ti, say all those 
corresponding to <3?(x) = P(x) e~ ax , where P(x) is an arbitrary polynomial 
and a is some fixed positive constant. Then denote by S* the space dual to 
S and, with the usual abuse of notation, consider Ti as a subspace of S*, so 
that we obtain the Gelf'and Triple (or Rigged Hilbert Space |23j) 



The space S* contains in particular the vectors |x') ® |±), where |x') is the 
generalized eigenvector of the Dirac x-coordinate operator q, 



corresponding in the coordinate representation to 5(x — x'). 

The introduction of the time interval At as in (fT^j) . in turn defines a 
length interval At on the x-axis (recall that c = 1 now) , and a corresponding 
direct-integral decomposition 



s < n < s* . 



(15) 



q\x') = x'\x') , 



(16) 




(17) 
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where V xo < S* is spanned by all vectors of the form \xq + mAt) (g) |±), with 
x G (—At/2, At/2] fixed, and m G Z. We note at once that each V xo is 
invariant under the action of the QRW evolution operator V U, i.e. 

VUV X0 <V X0 , (18) 

because 

V \x Q + mAt) <g> |±) = \x Q + (m =f l)At) ® |±) . (19) 

In order to describe the QRW evolution more fully, we now write the 
initial state with wave function as in (J5J) and (jZJ), as an entangled state of 
the walker and coin subsystems, 



^o» = / {c+(x)|x) <8) |+) + c_(x)|x) <g> |-)} dx, (20) 



where 

V27T J-oo 2y/E(p)(E(p) + 1) 

Normalization of |^o)) is satisfied because 

/oo 
{|c+(x)| 2 + |c_(x)| 2 } dx = l, (22) 
-oo 

as a consequence of (@J) and ©• At this point we emphasize again that 
although, as is well known [21], |x)®|±) is not a positive-energy (generalized) 
state, |^o) is a positive-energy state, as a consequence of the particular form 
of the coefficients in (|21|). 

The expansion in ()20j) can be rewritten as 

/At/2 
{c+(x + mAt) |x + mAt) O |+) (23) 
-At/2 



+c_(x + mAt) \x + mAt) <8> |— )} dx , 

which is to be compared with (|17p. In view of the invariance of each V XQ 
under the action of the QRW evolution, we can restrict our attention to that 
action on each substate 

l$*o)) = Yl i c +( x o + mAt ) l x o + mAt) (8) |+) (24) 

ra£ z 

+c_(s + mAt) \x + mAt) <g> |-)} v 7 ^ 
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with xq fixed, even though these substates are not normalizable, and are 
not positive-energy states. The point is that the general form of any such 
substate is preserved under the action of the QRW evolution V U, with no 
change in the value of xq. The inclusion of the multiplicative factor y/At in 
last equation is for later convenience with the normalization. 

Consider firstly the action of V on a general substate |$ )) G Vx , say one 
with Xq = for definiteness. We have 



V\$ )) = ^ |mAt)®{c + ((m-l)At)|+) + c_((m + l)At)|-)}v / At 



= ^2 ^[c Q (mAt) \(m + a) At) <g> \a)]y/At 

= (E + ® P+ + E„ ® P_) |$ », (25) 



where 



£±|mAt) = |(m± l)At) , P±|±) = |±) , P T |±)=0. (26) 

The action of U on |$o)) is easily seen from p2j) . which implies that 

U\+) = cos(At) |+) + sin(At) |-) , 

17 |-) = cos(At) |-) -sin(At) |+) . (27) 

Combining (|23J) and (J2*7|) . we see that 

V [/ |$ » = (#+ ® P+C/ + ® P-CO |$o» • (28) 

If we had taken f+(p) = 0, f-(p) = f(p) in (|7j), we would have written 
instead 

u (p) = 1 + E ^~P |+\ + j L+£M+£ |_\ (29) 

2 v /p( p )( j e;(p) + i) 1 2 v /p( p )( j e;(p) + i) 1 7 ' 

and we would have obtained 

/CO 
{c+(x)\x) ® |+) + c_(x)|x> ® |-)} ofrr, (30) 
•00 

where now 

c-{x) = -j=[ } + E ^+P f(p)e**dp. 
V^J 2 v /E(p)(E(p) + l) JKyJ 
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Then, decomposing |\&o)) into substates \& xo )) as before, we would have ob- 
tained on a state of this general form, say one with x = 0, that 

V C/|$o» = {E- ® P+U + E + ® P-U} |$ » • (31) 

We will treat here the first the second one can be treated similarly. 

To proceed we choose — n < < it, and set 

|0/At> = i- J2 e- im *\mAt) , (32) 



2tt 

m£ 



so that 



/7T 
e im *|0/At)d0. (33) 

Considering the evolution operator V U acting as in (jSH), but now with E± 
diagonalized, we have 

V U{<j>) = (e^P+ + e^P-) U . (34) 
The eigenvalues of this 2x2 matrix with parameter are 



A±(0) = cos cos At ± i a/ 1 — cos 2 cos 2 At . (35) 
Suppose that the corresponding eigenvectors are 

M0)) = / ++ (0)|+} + /+-(0)|-}, 

M0) = /_ + (0)|+> + /__(0)|->. (36) 

Then the eigenvectors of V U are of the form |0/At) <8> |f±(0)), with eigen- 
values A±(0). Expanding |$o)) i n terms of these eigenvectors of V U we get 



l$o» = / {^+(0) 10/ At) ® K(0))+s_(0) |0/ At) ® | v _(0)» VAt 

(37) 

where 

S±(0) = E (c + (mAt)/; + (0) + c_(mAt)/*_(0)}e^. (38) 

meZ 

Hence 

= (VC01*o»= r {^ + (0)A + (0f|0/At) ® K(0)) 

«/ — 7T 

+#_(0)A_(0) n |0/At)® |^_(0))}VAt d0. (39) 
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If we now denote by X n the random variable defining the "walker po- 
sition" after n evolution steps, then we obtain for the "quantum statistical 
moment" 

((X n ) k ) = (($ n \q k ® l|$ n )) = Tr s+T (|$n»P»|g* ® 1) 

= Tr T ((Tr 5 |$„,))(($ n |)g fc ) = Tr T (#g fe ) , (40) 

where the expectation value has been expressed in terms of traces over the 
translational degree of freedom (T) of the Dirac particle — the walker system 
in the parlance of QRW — and its spin (S) — the coin system for the QRW. 
This has allowed us to cast the "quantum statistical moment" in terms of 
the reduced density operator = (Tr s \§ n )}((®n\) which, as it provides all 
the necessary statistical information about the position of the Dirac particle, 
could have also been the main object of our mathematical investigation, as 
happens in most studies of QRWs. 

We proceed to determine the statistical moment of the translation oper- 
ator, which takes the form 



((Xn) k ) = / [gi{<P)K(m-^) k {9 + m+{<p) n } 

J — TT 

+^(0)Al(0)"(-^) fc {^-(0)A-(0) n }] ^(^) fe+1 > (41) 

or equivalent ly 



Hence asn-> oo, At — > 0, with t = nAt large, we have that 



+ |^(0)| 2 HA'_(0)A-(0)) ) gAf + OCnA*)*- 1 . (42) 



{{X n /nAt) k ) ~ r{|^ + (0)| 2 (-^(0)/A + (0)) fc (43) 



+ M0)f(-*X(0)/A_(0)) K j ^ At. 



This admits the following interpretation [25]: we can take as a ran- 
dom variable, a function Y from Q = S 1 x {+, — } to the reals, with Y = 
-iA' + ($)/A+($) on S 1 x {+}, and Y = -iA'_($)/A_($) on S 1 x {-}. Here 
$ : Q —>■ R is a random variable which projects on the circle S 1 with mea- 
sure \g + ((j))\ 2 At(d(j>/2ir) on S 1 x {+}, and measure |g_ (0)T At(d(j)/27r) on 
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S 1 x { — }. Since in the above limit, all the moments of X n /nAt agree with 
all the moments of Y, and the support of X n /nAt is compact, it follows that 
X n /nAt converges weakly to Y. Hence we have that 



lim P(yi < XjnAt < y 2 ) = P{y x <Y < y 2 ) 

nAt=t^oo 

= / t+ M0)I 2 a«| + / r M0)l 2 A t g, 



(44) 



where the intervals of integration are T± = y\ < (— *A±($)/A±(<&)) < y 2 . 
Then it follows that in order to determine the long time position distribution 
we need only determine g±(4>) and \±(4>). 

Suppose now that we specialize to the case of a highly localized initial 
electron state [H] with f(p) = f u (p) = e ~ p2 / 2l/2 / ( u^/ti) , where v is large 



and positive and quantifies the extent of the localization of the Dirac parti- 
cle's initial state — the larger is u, the sharper is the initial localization. As 
v approaches infinity one has that 



27T V \^ Jo 



c+(x) ~ -£= A /-^= / e wpx e- p2/2u2 dp, (45) 



o 



c(x) ~ -^=J-^= I e iupx e- p2/2,/2 d P . (46) 

V27T V V 71 " J-oo 

Note that in the limit u — > oo, f™. \c + (x)\ 2 dz = \c-(x)\ 2 dz = 1/2. If 
we now make uAt small by taking At small enough, then 

g ± (<j>) ~ ^2v^/(z/At2) e -* 2 /^ 2A ' 2 )/* + (0) i/ 0>O 

~ y/2y/Z/(vAt*) e^ 2 /^ 2At2 )/*_(0) i/ 0<O. (47) 

Also 

- *A' ± (0)/A ± (0) = ± /l Sin0C 2 °'f t 2 J = ( 48 ) 
V 1 — cos z Ar cos z 

Note also that |c/+(0)| 2 + (0) | 2 = 2 v ^Fe^> 2A '7(z/At 2 ). To compute the 
asymptotic distribution we need to compute the integrals 

27r ~ Ai F+W (2/)) 
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and 



,JK'\-m:-,A 27r 

1 /"^ 2 1 

- S Ejf l^^k-^t-,))! 2 ^^, (50) 

where the index % labels the local inverses of the function h. Because of 
()47|). the only inverse relevant to leading order is the one that keeps close 
to zero. For this inverse, h^ 1 (—y) = — (y). Furthermore, and realizing 
that \v±(fa))* = |f=p(— fa)), we can show that \g±(— fa)\ 2 = \g±(fa)\ 2 - A direct 
computation then gives 

At sin At r 1 hH P -(^(y)) 2 /^ dii 

2tt J yi (l-y^y/cos* At-?/ 2 ^At 2 

, (51) 

If we set y = h(fa), then for all inverses fa we have cos 2 0j = (cos 2 At — 
y 2 )/ cos 2 At(l — y 2 ). Since for our inverse the value of fa is small we have 



(V^)) 2 = 2 ~sin 2 ^ = l-(cos 2 At-y 2 )/(cos 2 At)(l-y 2 ) 

= ?/ 2 sin 2 At/cos 2 At(l - y 2 ) . (52) 

Taking the limit At — > we arrive at the asymptotic distribution function 
associated with the random variable X n /nAt ~ Y, 

P(Vi <Y<y 2 )= lim (h + I 2 )= / F(y) d, , 
v^/ti (1 — y i yi i 

In Fig. 1 we plot this two-horned probability distribution for three values of 
the localization parameter u, and recognize it as the 1-dimensional analogue 
of the result obtained for the Dirac particle in three dimensions in Eqn 
(3.1). 



4 Discussion 

It has been shown that the 1-dimensional Dirac evolution of a state with 
positive energy and definite spin is equivalent to a QRW, in the limit of 
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Figure 1: The asymptotic position probability density function, with local- 
ization parameter v — 1.9, 2.5 and 2.9. As v increases, the plots become 
more sharply peaked near the ends of the interval. 
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small positional steps and a large number of iterations. An initial state 
that is highly localized, with all but one momentum component set to zero, 
spreads in the remaining direction at a speed that almost surely approaches 
the speed of light as the initial localization increases. 

This relationship between the Dirac evolution and a QRW leads to the 
intriguing speculation that at some small space-time scale, there may really 
be a QRW defining the evolution of states of the relativistic electron, and 
that it is the Dirac evolution that is only a large scale approximation. One 
way to test this would be to make very precise measurements of the spreading 
characteristics of initially highly localized electron states over short distances. 
Comparison with the characteristics that are typical for a QRW, in particular 
the shape of the position probability distribution at early times, may reveal 
whether or not there is indeed a QRW underlying an approximate Dirac 
evolution. 

It is tempting to speculate further that there may be some deep rela- 
tionship between such an underlying QRW and the Zitterbewegung of the 
relativistic electron, as first discussed by Schrodinger This awaits fur- 
ther study. 
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Figure caption: 

Figure 1. The asymptotic position probability density function, with 
localization parameter v — 1.9, 2.5 and 2.9. As v increases, the plots become 
more sharply peaked near the ends of the interval. 
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